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Abstract. Let Q C K , N > 2 be a bounded smooth domain and a > 1. 
We are interested in the singular elliptic equation 

Ah= -h~ a -p inn 

a 

with Neumann boundary conditions. In this paper, we gave a complete de- 
scription of all continuous radially symmetric solutions. In particular, we con- 
structed nontrivial smooth solutions as well as rupture solutions. Here a con- 
tinuous solution is said to be a rupture solution if its zero set is nonempty. 
When N = 2 and a = 3, the equation has been used to model steady states of 
van der Waals force driven thin films of viscous fluids. We also considered the 
physical problem when total volume of the fluid is prescribed. 



1. Introduction 

The equation 

(1.1) h t = V • (h 3 v P ) 

has been used to model the dynamics of van der Waals force driven thin films 
of viscous 

fluids |2Z| EE] ES| EDI- Here h is the thickness of the thin film and the 

pressure 

(1.2) p=~h- 3 -Ah, 

is a sum of contributions from disjoining pressure due to attractive van der Waals 
force and a linearized curvature term corresponding to surface tension effects. 
Hence, becomes 

(1.3) h t = -V ■ (h-^h) ~V ■ (h 3 vAh) , 
which is a special case of the generalized thin film equation 

(1.4) ht = -V- (h m \7h) - V • (h n vAh) 

where the exponents m, n represent the powers in the destabilizing second-order 
and the stabilizing fourth-order diffusive terms, respectively. This class of equations 
occurs in connection with many physical models involving fluid interfaces|23 j [M | . 
For example, when n = 1 and m = 1, it describes a gravity driven Hele-Shaw 
cell[T]|3||ll | |14 | |2l) ] ; for n = m = 3 it describes fluid droplets hanging from a 
ceiling[i3|: and for n = and m = 1, it is a modified Kuramoto-Sivashinsky equa- 
tion which describes solidification of a hyper-cooled melt 0] [3- Over the past two 
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decades, these models have also been the focus of rigorous and extensive mathe- 
matical analysis^ |S1IE1D3]||I21P1||IS1I2D||2I1I221I2S|- 

As in the van der Waals force case, when n — to ^ 1 , letting 

(1.5) p= h m - n+1 - Ah, 

to — n + 1 

we can rewrite (II .4f) as 

fit = V ■ (h n Vp) . 

Now we consider viscous fluids in a cylindrical container whose bottom is rep- 
resented by 17, a bounded smooth domain in M 2 . Since there is no flux across the 
boundary, we have the Neumann boundary condition 

(1.6) ^ = on an. 

ov 

We also ignore the wetting or nonwetting effect, and assume that the fluid surface 
is perpendicular to the boundary of the container, i.e., 

dh 

(1.7) - = 0onSR 

dv 

Whenever m — n ^ —1 or —2, we can associate (11.4(1 with energy 

1 . — - .2 1 .ra-n+2 



E(h) = jf Q|v/i| a - 



(to — n + 1) (to — n + 2) 
and formally, using (|l.tj|) , (|1.7|) , we have 



m— n+1 1 



P v • (h n v P ) = - / h n \vp\ 



Hence, for a thin him fluid at rest, p has to be a constant, and h satisfies <|1.5ll . 
Therefore, letting a = — (to — n + 1), we are led to the elliptic problem 

Ah = \h~ a -p in Q, 
^ = Q a on dfl, 



(1.8) 



where f2 is a bounded smooth domain in Mr, N > 1 and p is a constant. 

When N = 1, this equation has been studied by R. Laugesen and M. Pugh in |2H] 
where they produced positive, smooth steady states for all a and touchdown steady 
states for a < 1. In [JJ], A. L. Bertozzi, G. Grim and T. P. Witelski considered ((1.8(1 
with additional Born repulsion term which leads to the elliptic equation 



(1.9) 



Ah=±h-<*(l-(lf)-p in n, 

§± = o on an, 



where (3 is a positive constant. When e > 0, the associated energy to (11.9(1 is 
bounded from below which makes a variational approach possible and enables them 
to show the existence of an energy minimizer in any dimensions. It seems difficult 
to extend this approach to the limiting case e = 0. 

The goal of this paper is to understand radial solutions of ((1.8(1 when N > 2 
and a > 1 which we will assume throughout this paper. In particular, when N = 2 
and a = 3, we come to the van der Waals force driven thin films in the physically 
realistic dimension. When a > 1, except the limited discussions in |18j . there 
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seems no established elliptic theory for (|1.8|1 . and hence it is the mathematically 
more interesting case. We remark that energy method can be applied to yield 
nontrivial solutions to (|1.8[1 when a < 1. On the other hand, the behavior of radial 
solutions is also quite different when a > 1. For example, we will show that for 
N > 2, the radial solutions will never vanish away from the origin which contrasts 
with the a < 1 case where touchdown steady states can be shown to exist in any 
dimensions. 

Due to the singular nature of l|1.8[) . we need to be careful in discussing "solu- 
tions" to Q1.8|l . We say h is a continuous solution of l|1.8f) in Q, if h ^ and is a 
nonnegative continuous function in fl satisfying the equation in Ijl.SJI in the open 
set {x £ f2 : h (x) > 0}. The rupture set of h, 

E = {ieS]:/!(i)=0}, 

corresponds to "dry spots" in the thin film, which is of great significance in the 
coatings industry where nonuniformities are very undesirable. Standard elliptic 
theory implies that h is smooth and hence a classical solution of (|1.8I) in 17\E. An 
interesting Hausdorff dimension estimate of £ can be found in |18j where it is shown 
that any finite energy solution satisfies W (£) = where /i = N — 2 + ^zj • For 
van der Waals force driven thin film, we have N — 2 and a = 3, hence Tt 1 (£) = 0, 
i.e., the thin film with finite energy can't have one dimensional rupture set. 
For any p > 0, let 

(1.10) £ = (op)"- 

then h = £ is always a solution to (|l.tj|) . The natural question is whether it is the 
only solution. For the radially symmetric case, after a simple scaling, the uniqueness 
theorem in |13| implies: 

Proposition 1.1. Let N > 2 and a > 1. For any given R > 0, there exists a 
constant pg, such that for any p < po, h = £ is the only radial solution of IjJ.iSjl in 
Br(0). 

When p is large, nontrivial solutions do exist. In fact, we have 

Theorem 1.2. Let N > 2 and a > 1. For any given R > 0, there exists a 
nondecreasing sequence of {pk}, such that for any p > pk, has at least k 

nontrivial smooth radial solutions in Br (0) . 

Theorem 11.21 is an application of Theorem 11.31 below which gives a complete 
description of all nontrivial smooth radial solutions. 

Theorem 1.3. Let N > 2 and a > 1. For any given p > 0, and for any rj > 0, 

T] there exists an increasing sequence {r^.' v }'^L 1 with 

where po is the constant in Proposition and 

(1.11) hm (^-rH^M^- 

smc/i £ftcrf for each r^' n ' , there exists a unique smooth radial solution of (1 i.ffl m 
-B r p.>? (0) satisfying h (0) = 77. 
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We say a continuous solution to 1)1.8(1 is a rupture solution if £ is not empty. It 
will be shown that for radial solutions, rupture can only occur at the origin. (See 
Corollary 12 . 21 below.) Our main result in this paper is as follows. 

Theorem 1.4. Let N > 2 and a > 1. For any given p > 0, there exists an 
increasing sequence | r fc'°| with 

(1.12) rf > 
and 

(1.13) ^(^-^^.(ap)-^, 

where po is the constant in Provosition M . 1\ such that for each r^'° ; there exists a 
unique radial rupture solution to in B rP ,o (0). Furthermore, if R ^ r^'° for 

any k, then there is no radial rupture solution to ^1.8t in Br (0). 

Remark 1.5. The rupture solutions constructed in Theorem \l-4\ above are weak 
solutions to l(i.#(l in the distributional sense. (See Remark \4-3\ below.) We also 
remark that when N = 1 and a > 1, there is no radial rupture solutions |2U| . 

For any given p > 0, there exists a nontrivial radial solutions to ((1.8(1 in Br (0) 
if and only if = i? holds for some 77 > 0, 7/ 7^ £ and for some integer fc > 1. 

In physical experiments, usually the total volume of the fluid is known, i.e., the 
average film thickness 

1 



h= w\L Hx)dx 

is given while the pressure p is an unknown constant. Hence, given h > 0, we need 
to find function h and constant p, such that 

( Ah = \h' a -p _ in Q, 
(1.14) I ^J n h(x)dx = h 

{ §1 = on an. 

When f2 = Si (0), all radial solutions of 1(1.14(1 can be obtained by scaling from 
solutions in Theorems 11.31 and 11.41 We will discuss such scaling in Section In 
particular, we will show 

Theorem 1.6. Let N > 2, a > 1 and 51 = B\ (0) C Mr. There exists a sequence 
of thickness hi, h,2, ■ ■ • satisfying 

lim Vknhk = 1, 

k — >oo 

such that for any k, I ] 1. 1$ with h = hk has a radial rupture solution which, viewed 
as a function in r, has exactly k—l critical points in (0, 1) . Furthermore, if h ^ hk 
for any k, then \1.14\ has no radial rupture solution. 

When £1 = B\ (0), Proposition 11.11 implies that nontrivial solutions to ((1.14(1 
must satisfy p > p^. Since 

P\ Jn a 

we may ask the existence of a critical average film thickness h$ so that there is no 
nontrivial solutions to ((1.14|) whenever h > h . Numerical analysis suggests that 
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ho does not exist. However, we are unable to provide an analytical proof. Such a 
proof could be possible if we have better understanding of h (p, 77, k) which is defined 
in section 5. Moreover, the detailed property of h(p,r],k) could also provide us a 
statement similar to Theorem II .61 for the smooth radial solutions. 

The paper is organized in the following way: In Section 2, we show that any 
radial solution can be extended to a global solution which is oscillating around 
£. In Section 3 and Section 4, we discuss smooth radial solutions and rupture 
solutions respectively and Theorems 11.31 and 1 1 . 41 are proved. In Section 5, we use 
scaling argument to prove Theorems 11.21 and 1 1 . (j| 



2. Preliminaries 

Recall that given a > 1 and p > 0, h £ C° {Br (0)) is said to be a continuous 
solution of 

(2.1) Ah = -h- a -p 

a 

in Br (0) if h > and it satisfies (|2.1(l in the open set {x £ Br (0) : h (x) > 0}. 

Let h be a radially symmetric solution to (|2.1|l . we can view h as a continuous 
function defined on [0, R) satisfying 

N - 1 

(2.2) h" + ti + f(h)=0 

r 

in the set S+ = {r e (0, R) :h{r)> 0}. Here 

f(h)=--h~ a +p 
a 

is monotone increasing and its antiderivative 

F(h)= 1 

a (a — 1) 

is convex in (0,oo). Let 

£ = (ap) " , 

then F' (£) = / (£) = 0, and F achieves its absolute minimum at £. Furthermore, 
F (h) — > 00 as h — ► 0^ or h — ► 00. 
For each r e 5 + , letting 

ei (r) = (r)) 2 + F(Mr)), 

e 2 (r) = i (r^ 1 /!' (r)) 2 + r 2 ^" 1 ^ (h (r)) = r 2 ^" 1 ^ (r) , 

we have 

(2.3) i.[ ei(r)] = _^l ( ^ (r )) 2 <o, 
dr r 

and 

(2.4) ^- [ ea (r)] = 2 (JV - 1) F (ft (r)) r 2Ar - 3 > 0. 

The monotonicity of e\ and e 2 will be used to obtain a priori bounds. 
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Lemma 2.1. Let < r± < r 2 be such that (ri,r 2 ) C S + . Given r 6 (ri,r2), rae 
have, for any r € (ri,r2) } 

(2.5) ei (r) < (-) ei (f) . 
Furthermore, 

(2.6) ci<fr(r)<c 2 

where ci,c 2 are two positive constants depending on a,p,N,rx,f,h(f) and hf (f), 
and are independent of r 2 ■ 

Proof. Since e 2 (r) is monotone increasing, we have, for any r € (ri,f], 

2(N-1) 

(2.7) ei (r) = r 2 ^e 2 (r) < r^-^e, (f) = Q e x (f) . 

On the other hand, since e\ (r) is monotone decreasing, we have, for any r € [f , r 2 ), 

(2.8) ej (r) < e x (f) . 

Combining (|2.7|) and (|2.8|) . we obtain (|2.5|) . Now for any re (ri,r 2 ), 

/ - \ 2(7V-1) 

F(/ l (r))<ei(r)< (-) e x (r) , 

so (|2.6|) follows from the fact that i 7, (/i) — > oo as h — > + or /i — > oo. □ 

Corollary 2.2. ft, can not have rupture away from the origin, i.e., S + = (0, R). 
Furthermore, h can be uniquely extended to a positive smooth solution of \2. 6 3, in 
(0,oo). 

Proof. Since S + is open, it is a union of open intervals of the form (ri,r 2 ) with 
fi,T2 4- S + . Given any such interval, if r\ > 0, Lemma 12 . II implies 

liminf h (r) > 0, 

and since h is continuous, we conclude h{r\) > 0, which contradicts the assumption 
r\ S + . Similarly, we can get a contradiction if r 2 < R. Hence, S + = (Q,R). 
Extending h to a maximal interval of existence (0, R*). If R* < oo, applying 
Lemma 12 . II again, we have for some positive constants c\, c 2 

ci < h (r) < c 2 for any r e (R/2, R*) , 

so the solution can be extended beyond R*. Hence, R* = oo. □ 

Now, redefining S + = {r > : h (r) > 0}, we observe that S + = (0, oo) and 
Lemma 12.11 still holds. In particular, l|2.6|l holds for all r% < r < oo. In the 
remaining part of this section, we shall show that h oscillates around £ near r = oo. 

We will need Sturm's Separation Theorem. 

Lemma 2.3. 17" Let q (t) be a real-valued continuous function such that 

< m < q (t) < M. 
Given ti > t\ > 0, if u = u (i) is a nontrivial solution of 

u" + q (t) u = 
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satisfying u (t) > on (t\,t2), then 



j4ncf z/ m addition u (t±) = u (£2) = 0, then 



M 



Lemma 2.4. For any ro > 0, there exists r± > ro smc/i i/ia£ h' (ri) = 0. 

Proof. Suppose this is false, then we have either h' (r) > for all r £ (ro,oo) or 
hi (r) < for all r € (ro, 00). Hence, h is strictly monotone increasing or decreasing 
on (ro,oo). From Lemma f2.ll and the observation above, it follows that h is also 
bounded at 00. So we can assume 

lim h(r) = ( 

r — >oo 

for some £ > 0. For any r > ro, integrating (J2.2|) from ro to r, we obtain 
(2.9) hi (r) - hi (r ) + 



N - 1 

/i'(s)ds + / f(h(s))ds = Q. 



Since 



\(h' (r)) 2 < ei (r)< ei (r ). 



hi (r) is bounded in [ro, 00). Now, as hi (r) does not change sign, 



hi (s) ds 



< 



N - 1 



ro 



hi (s) ds 



N - 1 



ro 



\h(r)-h(r )\ 



which is also bounded in [ro,oo). Hence identity (|2.9|l implies that 

f(h(s))ds 



is bounded in [ro, 00). Since linv^oo h (r) = £, we have 

lim f(h(r)) = f(0 

r— ►QO 

which must be 0. Thus £ = £ and 

lim /i (V) = £. 

r — >oo 

Now let 

N — 1 

i> (r) = r^~~ (h (r) — £) , 

then 

(2.10) v. 



— « + r 2 /(M = o. 



Since / (£) = 0, we can rewrite (|2.1U|) as 



(2.11) 
where 



B(r) = 



v rr + B (r) v = 0, 
f(h)-f(£) , (iV - 1) (3 - JN0 



4r 2 



Now linv^oo ft. (r) = £ implies 



lim S (r) = /' (0 > 0. 
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By Lemma \2. 31 

N — l 

v (r) = r~ ~ (ft (r) — £) 

will be oscillating around as r — ► oo, which contradicts the assumption that 
h(r) £ in a strictly monotonic manner. □ 

Next, we have 

Lemma 2.5. Let hi (ro) = for some r$ > 0. 

(i) . Ifh(r )=Z, thenh(r) = C 

(ii) . If h(r ) > then there exists T\ > r , such that hf (r) < on (r ,ri), 
h' (n) =0,h (n) andF (h (n)) < F (ft (r )). 

(raj. If < ft(r ) < £, then there exists r± > r , such that hi (r) > on (ro,ri), 
= 0, ft(n) > £ andF(/i(n)) < F(h(r )). 

Proof, (i). This is the standard ODE uniqueness result. 

(ii) . Since h (r ) > £, we have / (h (ro)) > 0. Now 

(r N -'h')' = -r*-if(h), 

implies that r*" 1 /!' is strictly monotone decreasing in (ro, ro + 6) for some <5 > 0. 
Hence we have h' (r) < on (ro, ro + 6). Applying Lemma ITU there exists r\ > ro, 
such that h 1 (r\) = 0, and we also have h! (r) < on (ro,ri) if we choose the 
smallest such r\. If h{r\) > £, we would have r N ~ x h' is strictly decreasing near 
r\, hence h! (ri) < 0, which gives a contradiction. And if h{r\) = £, then h = £, 
which contradicts the hypothesis h(r ) > £. Hence we have h{r{) < £. Finally, 
F {h (ri)) < F (ft (r )) follows from (IQl . 

(iii) . Similar to the proof of part (ii). □ 

Let ft- be a nontrivial global solution of (|2.2|) . starting with r\ > such that 
h' (ri) — 0. The existence of ri is guaranteed by Lemma 12.41 Without loss of 
generality, we assume h (ri) < £. For k = 1, 2, • • • , we define through Lemma 12.51 

r 2fc = sup{r > r 2 fc_i: h! (s) > for all s G (r 2fc _i,r)} , 
r 2 fc+i = sup{r > r 2fc : h! (s) < for all s £ (r 2fc , r)} . 

Lemma 2.6. 

lim r/c = oo. 

k — >oo 

Proof. If it is not true, then we have 

lim rfe = r* 

for some r* > 0. Since ft is smooth, we have 

&(»•*)=£, ft'(r*)=0, 
hence Lemma 12 . 51 implies ft = £, which is a contradiction. □ 

Next, we show that the lengths of oscillating intervals r k +\ — r k are bounded. 
Lemma 2.7. There exists positive constants C\ and C 2 such that 

C\ < {r k+ i - r k ) < C 2 

for any k = 1, 2, 3, • • • . 
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Proof. Since 



lim rfe = oo, 

k — >oo 



we only need to prove the lemma when is sufficiently large. Differentiating Ij2.2jl . 
we have 

h >» + ^Lzl h " _ Hzlh' + h- a - x b! = 0. 

r r z 

Let 

N — 1 

w (r) = r 2 h (r) , 

then w satisfies 

(2.12) w" + (a— 1 - «, = 0. 

Since h is bounded away from both zero and infinity when r — > oo by (|2.6|) , we have 
for some R > such that for any r > R, 

N 2 - 1 

d < fc— 1 (r) - < c 2 

for some positive constants ci, c 2 . Since r&, fc = 1,2, are zeros of u>, the 
conclusion follows from Lemma 12.31 □ 

Finally, we have 

Lemma 2.8. 

lim h (r) = £. 

r — > oo 

Proof. Starting with r*i > such that h' (ri) — 0, fo(ri) < £, we define as 
above. Since -F (h (r^)) is monotone decreasing in k, and /i (r 2 fe) > £, /i (^2fc-i) < £j 
the property of function i* 1 implies ft (r 2 fc) is monotone decreasing and h (r 2 k-i) is 
monotone increasing. Hence we have the limits 

r)i = lim h (r 2 fc) > £ > ^2 = lira ft (r 2 fc_i) . 

k — 'OO fc — >oo 

Now 

r r k+i 

Vi - m < \h (r fc+ i) - h (rfc)| = / (r)| dr 

<(r fe +i-r fe p \ti\ 2 d r y <Vc~ 2 ^ k+1 \h'\ 2 dr 

which implies 

(2 . 13) r' M! * > _L P+' M . > (*-■»)' > (*-»>' r- ^ 



r 



when fc is sufficiently large. In the last inequality we used rk+i < 2r^ when fc is 
large. From ()2.3|) . we have for any r > ri, 

[ei (ri) - ei (r)] < 



r iV — 1 WJ ~ iV - 1 

Therefore J— L- is integrable at oo. Since - is not integrable at oo, (|2.13|) implies 
that 771 = 772 = S, and 

lim h (r) = £. 

r — >oo 

□ 
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Corollary 2.9. 

(2.14) lim (rjfe+i - r k ) = tt (ap)~^ . 

k — >oc 

Proof. In equation l|2.12l) . we now have 

/ n AT 2 -l\ , i+c. 

lim h^" 1 - = e— 1 = op)" . 

r—>oo \ 4r / 

Hence (|2.14() follows from Lemma l2~3l □ 

3. Nontrivial Smooth Radially Symmetric Solutions 

This section is devoted to the proof of Theorem 1 1.31 
Given 77 > 0, we consider (|2.2|l with the initial values 

h(Q)=r)> 0,h r (0) = 0. 

The local existence and uniqueness of a smooth solution is standard since / is 
smooth when h is bounded away from zero. And such solution is actually a global 
solution from Corollary 12.21 For any 77 ^ £, without loss of generality, we assume 
77 > £. Since 



(r N - l ti)' = ~r N - l f{h) = -r"- 1 \-^h 

we have (r^ -1 /),')' < in (0, S) for some small S > 0. From h! (0) = 0, we conclude 
h' (r) < in (0, 8). Then we can define 

7*1 = min {r > : b! (r) = 0} . 

The existence of T\ > is guaranteed by Lemma l2~^l with ro = | . From the analysis 
in the previous section, /i will be oscillating around £, and all critical points of h 
can be listed as r± < r-i < < ■ ■ ■ , with 

C\ < r k+ i - r k < C 2 , 

and 

lim (rk+i - r fc ) = 7r (ap) 2a . 

k — >oo 

Hence for any k > 1, h is a nontrivial smooth solution of 

Ah=±h- a -p in B Pfc (0), 
f = on 5B rfc (0). 

And all nontrivial smooth radial solutions of l|1.14|) . when £1 is a ball, can be 
obtained this way. More precisely, let f2 = B R (0) for a given i? > 0, then (|1.14|) 
has a nontrivial smooth radial solution if and only if R — for some 77 > 0, 77 7^ £ 
and for some k > 1, here we write r/c = r^' 1 ' to recognize its dependence on p and 

77. 

Now we recall the uniqueness result of M. Del Pino and G. Hernandez |13| . 
Proposition 3.1. Given a > 1, there exists oIq > 0, such that 

-d/\u + v,- a = l in Bi(0), 
§f=0 on dBi{Q) 

has no nontrivial radial solution whenever d > do- 



1 



.-a 



(3.1) 
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It is easy to verify 

h (x) — (ap)" h (rix) 

satisfies (|3.1[) with 



d=- " 



(ap) " r\ 



Hence Proposition 13 . 1 1 implies 
ri> ( 

here 



a + l 

,(ap) " d , 



+1 

Po ' 



po = ( ) = a *+id a+1 . 
\ a ° d J 

In general, depends on both p and r\. We refer to Corollarv l5.2l for the scaling 
of rj; when p, r\ changes. 

Lemma 3.2. For any r\ > 0, 77 7^ we /lave 
In particular, 



2Na (r) - £) 



lim ri (77) = 00. 

r\ — >oo 

Proof. First we assume rj > £. From the definition of r\ (77) and Lemma 12.51 we 
have h! (ri) = 0, h (ri) < £ and for any r e (0, ri), < ft, (r) < 77, /i' (r) < 0. Now 

/ t—OL h~ a \ £ — a — a 

(r^/V)' = ( S — ) > -r^ 1 - — ■ 



a a 



Integrating from to r, we have 



r a (r) > 



AT a 
i.e. 



fc'(r)>~(r--0 



Integrating again from to r\ , we have 



Mri) - MO) > (r a -»T Q ) 



hence 



^ (??)> / 2iVa (ft (0) - h (n)) > / 2iVa (77 - 



£ — a T}~ a — y £ _ a — Tj~ a 

The bound when 77 < £ can be proved similarly. □ 



12 



HUIQIANG JIANG AND WEI-MING NI 



4. Rupture Solutions 

In this section, we will consider radial solutions to H2.1(l which are not smooth 
and prove Theorem 11.41 From Corollary 12.21 we need to consider h £ C° ([0, oo)) 
such that h (0) = and h satisfies 112.2ft in (0, oo). 

First, we check the growth rate of h near the origin. 

Lemma 4.1. Let h be a radially symmetric rupture solution, then for any S > 0, 
there exists positive constant c\ such that 

h (r) > c\r "+ 1 

holds for any r € [0, 5]. 

Proof. Since h is positive and smooth away from the origin, we only need to prove 
the bound for small 6. Let S > be sufficiently small so that 

-h~ a (r)-p>^-h- a (r) 
a la. 

holds for any r G (0, 6]. Now 

(4.1) (r^h'Y^r^^-p^^h-^ 

implies r w_1 /i' is monotone increasing in (0, 6]. Since h (0) = and h (r) is positive 
away from the origin, there exists a sequence j-j — > such that hi (r,) > 0. Hence, 
h' (r) > for any r € (0,5]. Integrating ()4.1|) from e to r, and using the fact that 
h is increasing, we have 

r*-i h > ( r ) - e N ^h> (e) > ^h^ (r) (r N - e N ) . 



Letting e — > 0, we have 



Hence for any r £ (0,5], 



r^h'M > —h- a (r)r N . 
V ; - 2Na y ' 



d . a+l ^ a + 1 



— (r) > 



Integrating from to r, we have 



i.e. for any r S (0, 6], 



h(r)> r°+ 1 . 



□ 



Lemma 4.2. Lei /i be a radially symmetric rupture solution, then we have for some 
positive constant Ci, 

h (r) < C2r a + 1 

for any r £ [0, oo). 
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Proof. Since h is uniformly bounded at oo, we only need to prove the inequality 
near the origin. First we claim 

(4.2) lim /^/i' (r) = 0. 

J-->0+ 



From 



it follows that r N 1 h! is monotone increasing near the origin. Thus, if l|4.2[l is false, 
we would have 

r N ~ l h'(r) > c> 

near the origin, hence 

ti (r) > cr 1 ~ N . 

Since r 1_Ar is not integrable near zero, the above inequality contradicts the fact 
that h is continuous. 

Given 5 > 0, for any r G (0, 5), 

JV-IlA' _ „AT-1 / " „ \ ^ n „iV-l ^ c l „N-l- 



(r"- L h'Y = r N - L p)< — r N - 1 < ^ r ^-i-^r 

\ a J a a 

by Lemma 14.11 Integrating from e to r, we obtain 

r"-^' (r) - e""^' (e) < ——^ -c^ a (r N ~^ _ 

«(^-^r) V 

Letting e — > 0, we have 

r^/i' (r) < — r cr a r w -^ 

a(JV- 

\ a+l 

i.e., 

ft' W < r cr Q r 1 -^r. 

Integrating from to r, we have, for any r S (0, (5), 

~l~ 1 _ 2 

h (r) < - r-c 1 Q r°+i . 

2« 



□ 



Lemmas 14 . 1 1 and 14 . 21 imply that (r) is of order r Q+1 near the origin. Now we 
write 

2 

(4.3) /i = cV(r)r° +1 , 

where 

2a ( 2 
N-2- 



q + 1 \ a + l 

Observe that h — c*r^T^ is a solution of 

Ah- -h~ a = 
a 
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in (0,oo). Direct calculation yields 

(4.4) ^ + {A + 1) ^_ + 9_M + Cr -^ =Q 
where 

A = iV-2+— ^- >0, C= — >0 
a + 1 c* 

and 

(?) = -H- (jv - 2 + -i-) - <p-°) . 

a + 1 \ a + 1 J 

Lemmas 14 . II and 14 . 21 imply that 

(4.5) < liminf ip(r) < limsupy(r) < 00. 

r^0+ r- 1-0+ 

On the other hand, let ip be a positive solution of 1)4.4(1 satisfying 1)4.5(1 . Then h 
defined by 1)4.3(1 is a rupture solution. 

Locally, there exists at least one solution of 1)4.4(1 with initial values 

(4.6) <p(0) = l,<pf(0)=Q. 
To see this, we rewrite the equation as 

/+( , +1) ^ + m^ + £Mz4)M +Cr -* = o 

Denoting 

%l) = tp-\, 

we have 

^ + (A+i)^ + ^ + £M + c7r -^ = o, 

where 

satisfies g (0) = 0, g' (0) = 0. Now let 01, a 2 be two numbers satisfying 

ai +a 2 =A = N-2 + ^—, a 1 a 2 =g / (l) = 2(N-2 + —!- 

a + 1 \ a + 1/ 

then the real parts of a\ , a 2 are both positive and it is easy to verify that 
Hence, we have 







Let 

Lip = —. ^7 r^-r~ ai I Is^-^-'l t a ^ 1 g (ijj (*)) <ft [> ds , 

( a i + ^r)(« 2 + ^r) y ° 1 y ° 

then for <5 sufficiently small, L is a contraction mapping from 

X = {V> S C([0,<5]) : |V»(r)| < 5 for any r G [0,(5]} 
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into itself. Here L is a real mapping even though oi, 0,2 could be complex numbers. 
Let ip be the unique fixed point of L in X, then ip = 1 + xj) is a solution to (|4.4(l 
satisfying Q4.6|) . 

Let t/3 be the local solution of (|4.4() we just constructed, then ft, defined by l|4.3[) is 
continuous with h (0) = and satisfies (|2.2() in (0, 5). Such solution can be uniquely 
extended to a solution in (0, 00) which converges to £ by Lemma l2.8l Thus we have 
constructed a global rupture solution. 

Remark 4.3. From the bounds in Lemmas \4-1\ and \4-S[ it is easy to see that the 
rupture solution we constructed is actually a weak solution of \2. 1\ in 10 x 
More precisely, we have 



and 



h € W& (M ) , ft~ Q S Lf oc (M ) /or any 1 < P < 
ftA0 = / f iftr a -p) 



ftoZds /or any G Cq° (R*) . 

From the proof of Lemma f4. II we have ft' > near the origin, so we can define 

ri = min {r > 0: ft' (r) = 0} , 

the existence of r\ is guaranteed by Lemma [2.41 Furthermore, as in the smooth 
solution case, we have a sequence {ffc} fe=1 such that for each k, ft is a rupture 
solution of 

Ah=±h~ a -p in B rk (0), 
f = on 0£ rfc (O). 

In the remaining part of this section, we will show that the rupture solution to 
(|2.2() is actually unique. 
In 14. 4j) . with 

r = e~, cj){t) = ip{r), 

direct calculation yields 

(4.7) tt - ^L0 t + 5 (0) + Ce^' =0 
on (—00, 00). 

Lemma 4.4. There exists a unique global solution to \4- 1\ satisfying 

(4.8) < lim inf (t) < lim sup (£) < 00. 

t— >°° t— >oo 

Noticing that r — > + is equivalent to i — > 00, the uniqueness of rupture solution 
follows from Lemma 14.41 Before proving Lemma 14.41 we first study the behavior of 
at 00. 

We write 

A2 AX- 



G (0) = [N - 2 + , , 

KV> a + l\ a + l J \2 a- I 

hence, 

G' (0) = <?(0). 
Multiplying equation l|4.7|l with t , we have 

(4-9) id + G (0)") = Atf - Ce"** > ±£ - f-e 



dt\ 2 v 7 1 ~ 2 ™ 2A 
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Lemma 4.5. Let cf> be a global solution to | |^.7| | satisfying then 
— oo < liminf t (£) < limsup^t (i) < oo 



Furthermore, the limit 



exists and is finite. 



b\ < oo. 



lim I & +G(d>) 

t — >oo \ Z 



Proof. If 4>t is unbounded at oo, then +G (<p) will be unbounded, so there exists 
a sequence {tk}kLi with lim^oo = oo, such that 



k — >oo 

For any t > tk, integrating 



lim >^M + G{(j)i tk))) = x. 



from ifc to i, we have 



(*) , tf(t k ) 



+ G (0(t))> 



2 



+ G(0(t k ))- C^e-^ds 



>—- + G(cp(t k ))-— ^— . 

Hence, 

lim(f + G W )=oc. 
From Q4.8|) . G (0) is bounded at oo, so we deduce 

hm — = oo, 

t — *oo 2 

which is impossible for bounded cf>. The L? (0,oo) bound of <f> t follows from (|4.9|) 

Ay 2 

and the fact that % + G (</>) is bounded at oo. Finally, since the right hand side of 
is absolutely integrable at oo, we have for any to, 



lim 

t — >oc 



,2 



0? (*o) 



+ (*„))+/ U^-Ge 



A2 ^„ S+T 



which is finite. □ 
Lemma 4.6. // 

lim (/) = <;, 

t — >oo 

then 

<T = 1. 
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Proof. If <r ^ 1, we first assume ? < 1, then for some small S > 0, 

for any t > To, where To is a sufficiently large constant. Hence 
(4.10) (e~ M 4>t) t > Se~ At 

for any t > T . Now since <pf is integrable, we can choose T\ > T with |0 t (7i)| 
sufficiently small. For any t>T\, integrating (|4.1U|) from 7\ to t, we have 



Hence 

<k (*) > 



2A 



when i is sufficiently large, which contradicts the boundedness of <f> t at oo. The 
case <; > 1 can be treated in the same manner. □ 



Lemma 4.7. 

and hence 



lim G {<!>{£)) = G(1). 

i — »oo 



lim </>(t) = 1. 



Proof. Since G (1) = min G (0), if 



lim [f + gm j : Gil). 



then lim t ^oo G (0) exists and equals G (1), as desired. We proceed by contradiction 
and assume that 



We claim 
Otherwise, 
which implies 



lim [4+G(^)=L>G(l). 



liminf G{4>) < L. 

t—*oo 



lim G {(f) 

t — >oo 



lim 4> = <; 

t — >oo 



for some <^ with G (<,) — L, a contradiction to Lemma [4.61 Hence, there exists a 
sequence {tkW^Li such that — > oo and 

G(l)<G(cf>(t k ))<L-6 

for some 8 > 0. Now we consider 

s fc = sup js > i fc : For any t <G (tfc, s) , G (0 (t)) < L - ^ 
Observe that s k is finite. Otherwise 0j > | for any t sufficiently large, and then 



is monotone with derivative bounded away from zero, hence it will be unbounded, 

8 

2 - 

\tf>(t k )-<f>(a k )\>5 1 , 



which gives a contradiction. Since G (cj) (s k )) = L — |, we must have 



18 



HUIQIANG JIANG AND WEI-MING NI 



where Si > is a constant depending on 5, L and G. Since 4>t is bounded, we have 



\s k -t k \>5 2 = -^- >0. 



However, for t% so large that 



f + G^>L- S - 



on (tfe,oo), we have 
on (t k , Sfe). Hence 

On the other hand, Lemma says 



0? W > | 



> \^ 



lim / 2 =O, 

which is a contradiction. □ 

Now we are ready to prove Lemma 14.41 

Proof of Lemma \J^\ Let <f> and <f> be two global solutions of (|4.7(1 satisfying 114.811 . 
Letting ip — <p — <f>, we have 



Here 



and 



Since 



we have 



ip u - Aipt + B (t) V = 0. 

4 



A = N- 2 



fl (t) - fiV - 2 ■ 2 



a+1' 



-a A~cn 



a + 1 V a + 1 



lim = lim (f>(t) = 1, 



lim S (t) = B = 2 ( AT - 2 + — ) > 0. 



t^roc \ a + 1 

It is easy to check that for any A such that 

A 2 - AX + B a = 0, 

we have Re A > 0. Since ip is bounded at oo, Lemma 14 . 81 below with Ao = implies 

1/3 = 0. □ 

The following result seems standard and should be well-known. A proof is in- 
cluded here for the convenience of the reader. 
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Lemma 4.8. Let u satisfy the linear equation 

(4.11) u tt ~ Aut + B (t) u = 0. 

Here A is a constant, and B (t) is a continuous function such that 

lim B(t) = B Q . 

t — >oc 

Let Ai, A2 be solutions of 

A 2 - AX + B a = 0. 
Suppose that there exists a constant Ao satisfying 

Ao < A m = min (Re Ax, Re A2) 
such that, for some positive constants T and c, 

\u(t)\< ce Xot 

holds for any t >T, then u = 0. 

Proof. Let u be any function satisfying (|4.11[l . For any A £ (Ao, A m ), let v 
It is easy to check 

(4.12) v tt -(A- 2A) v t + (B (t) -AX + A 2 ) v = 0. 
Since A < X m , we have 

A - 2A = Re Ai + Re A 2 - 2A > 2A m - 2A > 

and 

B a - AX + A 2 > 0. 
Multiplying (|4.12|l with vt, we obtain 

j t (vf + (B -AX + A 2 ) v 2 ) =2 (A- 2A) v 2 + 2(B -B (t)) vv t . 
For any s\ > since 

there exists T\ > 0, such that 



lim B(t) = B , 

t — ^OO 



\{Bq-B (t)) ««t I < ei (w 2 + (So - AX + A 2 ) w 2 ) 
holds for any t > Ti . Hence for any t > Ti , we have 

j t (v 2 + (B -AX + A 2 ) v 2 ) <2(A-2X + £l ) (v t 2 + (B -AX + A 2 ) 
Gronwall's inequality then implies that for any t > Ti, 
(4.13) Wt 2 + (B - AA + A 2 ) v 2 < c^e 2 ^" 2 ^)* 

where 

c £i = [ v 2 (Ti) + (Bo -AX + A 2 ) « 2 (TO] e - 2 (^- 2A +^) T 1 . 
Now let ui be the solution of (|4.11l) in the Lemma such that 

\ui {t)\ < ce Xot 
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holds for any t > T. Then for any A £ (Ao,A m ) and for any e £ (0,A— 2A), 
Vi = e~ xt u\ satisfies 

| (vl t + (Bo -AX + A 2 ) vj) 

=2 {A - 2A) vl t + 2{B -B (t)) v lVl}t 
>2{A-2X-s) vl t -vf 

for any t > T £ if we choose T £ >T sufficiently large. Since A — 2A — e > and 

(4.14) <ce- (A - Ao)t 

holds for any t > T, we have, by similar arguments as in the proof of Lemma l4.5l 

\vi,t (t)\ <Ci 

holds for some positive constant Ci and for any t > T e . Hence for any A S (Ao, A m ), 

|u!, t (t)| = |e At « 1:t (t) + Ae At Wl [t)\ < C 2 e At 

holds for some positive constant C2 and for any t > T e . Especially, for any 62 £ 
(0, A m — Ao), we have 

|ui,t {t)\<C 2 e^'+^ t 

holds for any t > T e where C2 is a large constant. Now for fixed A £ (Ao, A m ), we 
have 

(4.15) Kt| = |e- At u M - Ae- At ui| < C 3 e {x °- x+£2)t 

for any t > T e where C3 is a large constant. If u\ is a nontrivia! solution to (|4. 1 1|> . 
then vi is a nontrivial solution to (|4.12|l . Let V2 be another solution of (|4.12|l which 
is linearly independent of v\ . Then (|4.13|) holds for v 2 ■ Combining with (|4.14|) and 
l|4.15p . we have for any t > max{Ti,T e }, 

W (t) = VlV2,t - 

satisfies 

(4.16) \W(t)\<Kt+^) i Kt+4) i 

<( ^< 4e (Ao-A+e2)t . e (A-2A+£i)t _ q JA-2X+X - A+ei +e 2 )t 

On the other hand, since 

W (t) = {A- 2A) W (t) , 

we have 

(4.17) W (t) = W (0) e {A - 2X)t 
Choosing e\ and £2 small enough so that 

A - A + £1 + e 2 < 0, 

we conclude from 1)4.16(1 and ()4.17|l that W (t) = which contradicts to the as- 
sumption that Vi, V2 are two linearly independent solutions. Hence vi = and 
ut = 0. □ 
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5. Scaling of solutions 

In this section, we will use a scaling argument to prove Theorems O and Ol 

— — 

Let h p ' v be the unique solution to l|2.2|l satisfying h (0) = 77 ^ (op) " . When 
i] = 0, h p, ° is the unique rupture solution. Let r?'' J , k = 1, 2, • • • , be the increasing 
sequence of positive critical points of h p ' v . Then 

h*** (x) = (rf 7 )" 1 ^ h™ {t*'"\x\) 

satisfies 

f Ah = \hr a -pP^' k in Si(0), 
( f = 0° on 9^(0) 

with 
Let 

ft (p, V, k) = —\- I h*"'" (x) dx = f h*" (x) dx. 

\ B l\V)\ JB 1 {Q) \ B r k JB rh (0) 

Then h p,v,k (x) is a solution to l|1.14[) with 

h = h(p,rj,k) . 

_ j_ 

Lemma 5.1. For any p > and 77 > 0, 77 ^ (ap) a , 

ft"' 1 ' (a?) - (ap) 4 ft* 'to)*" ({ ap )^ x\ . 

Proof. Let 

/ ( x ) = (ap)~« /7^^ q p)"" ((<*p) ^ ^ , 
we have / (0) = 77 and 

A/ (z) = ap f A/i^'^M ({ap) 1 ^ a) 

I a \ / a I a 

So the lemma follows from the uniqueness of the radial solution. □ 
Corollary 5.2. For each k, 

and 

h(p,r),k) = h (ap)~ -q,k^ . 

Now we are ready to prove Theorem II. 21 

Proof of Theorem M.SX Lemma |2.5I implies that for fixed 77 > 1, r£' , = 1, 2, • • • , 
are well seperated. Hence, we can apply standard ODE theory to conclude that 

— v 

for each k = 1,2, rg' , viewed as a function of 77, is continuous in (l,oo). 
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From Lemma |3.2I and noticing r?' is monotone increasing in fc, we have for each 



fc = l,2,- 



lim rff' v 

r] — >oo 



Rk = inf 



Let 

v>i 

then Rk is monotone nondecreasing in fc. Furthermore, the interval (Rk,oo) is 
contained in the range of r£ ' . Given R > 0, let 

1 fRk\^ 

then pk is monotone nondecreasing in fc. For any p > pk, and for any 1 < i < fc, let 
r/i G (i?i,oo) be such that 

Then we have 

r *,(<xp)-im = {ap) -^ r hv i=Rf 

i.e., ft, p '( ap ) = is a nontrivial smooth radial solution to (|1.8|) in £?^(0). Since 
h vAcp)-«Vi : viewed as a function in r, has exactly i — 1 critical points in (0, R), we 
have found fc distinctive solutions as desired. □ 

- v 

Remark 5.3. For fixed k > 1, numerical computation suggests that rS' is not 
monotone increasing for 7/ £ [0, 1) U (1, oo). Hence, given p > po, we may have two 
different solutions with the same number of critical points in (0, R). 

From Corollary 15. 21 to get a solution to l|1.14|) through scaling, we can fix either 
p or r/. Without loss of generality, we assume p = — , and 77 1, 77 > 0. 

Proof of Theorem M.bl All radial solutions can be obtained by scaling. So (|1.14|) 
has a rupture solution only when 

h = hk = h ( — , 0, fc 

Let 
then 



i 



(r*r* 



(x)da;. 

B,.(0) 



Since 
we have 



|5 rfc (0)| 

lim (r) = 1, 

r — >oo 

lim -j / ft"' (2) = 1. 

= (0) 



|s rt (o)i y Br 

Hence the conclusion that 

lim \/k~Trhk — 1 

k — >oo 

follows from Corollary 12.91 □ 
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